Global Level Number Variance in Integrable Systems 
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We study previously un-researched second order statistics - correlation function of spectral staircase and 
global level number variance - in generic integrable systems with no extra degeneracies. We show that the 
global level number variance oscillates persistently around the saturation spectral rigidity. Unlike other second 
order statistics - including correlation function of spectral staircase - which are calculated over energy scales 
much smaller than the running spectral energy, these oscillations cannot be explained within the diagonal ap- 
proximation framework of the periodic orbit theory. We give detailed numerical illustration of our results using 
four integrable systems: rectangular billiard, modified Kepler problem, circular billiard and elliptic billiard. 



I. INTRODUCTION 

Interest in semiclassical properties of classically integrable 
systems picked up recently with deeper understanding of large 
persistent oscillations of the level number variance over an en- 
ergy interval as a function of the interval width and of the phe- 
nomenon of level repulsion as manifested through deviations 
from the Poisson statistics of nearest level spacings. llj-|5|] The 
precise nature of these effects are revealed by the structure of 
the correlation function of the level density and all relevant 
quantities can be computed both from the periodic orbit (PO) 
theory @| and by direct quantum-mechanical calculation 0]. 

An early attempt of evaluation of the correlation function of 
spectral staircase (SS) I7|] and interest in global level number 
variance (GV) were motivated by the fluctuations of thermo- 
dynamic quantities of mesoscopic electronic systems [8[]. For 
instance, evaluation of orbital magnetic response in the inte- 
grable circumstance |9l,[l0|] calls for performingensemble av- 
eraging (achieved via parametric averaging QjJ, |2j |4[) prior to 
thermal averaging, which, in turn, requires knowledge of the 
magnetic field dependence of the correlation function of SS. 

The central result of this work is to establish, theoretically 
and numerically, that GV exhibits large persistent oscillations 
around the saturation spectral rigidity. [[JJ 0-eI Moreover, 
it is shown that these oscillations cannot be described in the 
standard framework that utilizes the diagonal approximation 
(DA) of the PO theory |6|] but rather require an account of 
interference between periodic orbits with different winding 
numbers. Additionally, we evaluate the correlation function 
of SS and show that it can be expressed in terms of interval 
level number variance. 

This paper is organized as follows. In Sec. II, we evaluate 
GV and the correlation function of SS using the PO theory. 
In Sec. Ill, we present numerical evaluation of the GV vis-a- 
vis the saturation spectral rigidity for the rectangular, circular 
and elliptic billiards (RB, CB, EB) and the modified Kepler 
problem (MK). For RB, we proceed with a more extended 
analysis of the SS, its correlation function and interference 
effects in GV. 



II. THEORY 

A. Correlation Function of Spectral Staircase 

In PO theory, the fluctuating part of SS, 8Jr"{e) = jY{e) — 
,yV(e)), is found as a sum over POs and their time-reversals 
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where fi = (N — 1 ) /2 with N the dimensionality of the po- 
sition space |6| and 8j = 1/2 if the PO and its time-reversal 
coincide and 1 otherwise. Fluctuations of SS and fluctuations 
of level density 8p (e) are related via |6|] 
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with the use of 7) = dSj/de and on the account of the fact 
that the dominant contribution comes from differentiation of 
the oscillating term. 

The correlation function of SS is found from (£TJ as 

K^(e,co) = (8^(£ l )8^(£ 2 )) 
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where e = (ej + £2)/2 and CO = £2 — £\ C£. For integrable 
systems, ensemble averaging is understood as the parametric 
averaging |[l], 0, |4J] and the second, rapidly oscillating cosine 
was dropped in ([3]) as it produces a negligible contribution 
upon such averaging insofar as ©-dependence is concerned. 

In what follows, unless explicitly stated otherwise, we drop 
the argument of Aj and 7). We notice that the (interval) level 
number variance over the energy interval of width co is given 
byli 
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Consequently, combining (0 and (0, we have Q 
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where A is the mean level spacing and A"(e) is the saturation 
spectral rigidity given by 131 
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Below, it is ordinarily assumed that CO > 0. As expected, dif- 
ferentiation on £12 inside the cosine of the last equation of 
(0, also gives the correlation function of the level density 121, 
which can also be obtained directly from (0. 



B. Global Level Number Variance 



We now turn to GV of SS, which is defined as follows: Q 
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Formally, GV is a particular case of K^y(e,Co) in (0 with 
CO = 0. We point out, however, that (0 was obtained using 
DA, which is sufficient for evaluation of the ft)-dependence of 
Kyje.co) for a given e (see supporting numerical evidence 
below). IU2I1 However, it breaks down when the e -dependence 
of E ? (e) is considered. Whereas DA yields, upon averaging, 
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interference between POs in accordance with (0 
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must be considered to account for the full e-dependence of 
Eg(e) obtained in the numerical calculation below - namely, 
the persistent oscillations of E s (g) around A~(e). 

Indeed, the off-diagonal contribution to GV contains, per 
(0, the following term: 
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where the dependence of the action on the system parameter 
a is explicitly indicated. Parametric averaging involves inte- 
gration over the distribution function p(oc), such as a Gaus- 
sian distribution centered around the central value (Xq. ]4|] 
Ordinarily, such integration with rapidly oscillating terms in 
([Tol l produces negligible contributions (as is the case with the 
dropped term in (0 and (0 as well as with the off-diagonal 
contribution when CO ^ 0). The notable exception occurs when 



d a (Sj(s,a) ±Si(e,a))\ a =a Q = and the arguments of the 
cosines scale as (a — (Xq) 2 . Below we illustrate this circum- 
stance on RB. 

1. Rectangular Billiard 

For a particle of mass m in a RB with sides a and b, the 
amplitude, period and action of a PO with winding numbers 
M = (Mi ,Mi) are given respectively by H 
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Consider, for instance, interference terms between M = 
(M U M 2 ) and M p = (M 2 ,M\). Setting fi = 1 for simplicity, 
the cosine arguments in (fTOb have the following form: |Q] 0] 
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where a = a 2 /b 2 is the aspect ratio of RB. Ja] It is then triv- 
ially seen that at Oo = 1 and for Mi ^ M 2 , d a S- ^ while 
d a S+ = 0. 

Parametric averaging is performed via integration with the 
Gaussian distribution function p (a) whose width is <C 1, cen- 
tered at Oo = 1 . [3] Clearly, for a square, M and M p represent 
the same orbit per a 90° rotation; for RB with aspect ratios 
close to unity, that is a near square shape, we observe inter- 
ference from geometrically similar orbits with nearly equal 
lengths. We emphasize that since it is the near equality of 
lengths that matters, this argument can be easily extended to 
an arbitrary aspect ratio. 



III. NUMERICAL RESULTS 

In what follows, we express all energies in units of the mean 
level spacing A by setting A = 1 . 

A. Global Level Number Variance 



In Fig. we plot L g (e) vis-a-vis A"(e) UM for RB, MK, 
CB and EB respectively. IU4ll We observe that A~(e) ~ ^/e 
in RB. HI la] In MK, the saturation spectral rigidity exhibits 
quantum jumps to higher plateaus, while it experiences an 
overall growth as A~(g) - e 1 / 3 . Jlll3] In CB and EB, while 
scaling overall as A"(e) ~ y/s, as expected in a hard-wall bil- 
liard, the rigidity exhibits a far more complex behavior than 
in RB. 1 5] While not fully understood, we speculate that its 
origin may lie in the coherent effects of type-R orbits Hllll 
of approximately equal length - or length multiples - giving 
rise to global fluctuations of the level density. ^ IU5I1 As 
was already mentioned above, persistent oscillations ofL g (e) 
around A3 observed in Fig. cannot be explained in the DA 
framework. Below, we concentrate on RB in order to explicate 
the nature of these oscillations. 



Rectangular billiard 



Modified Kepler 




FIG. 1: Comparison between global level number variance and saturation spectral rigidity. Blue line: L g (e) of rectangular 
billiard, modified Kepler problem, quarter circular billiard, and quarter elliptic billiard. Black line: A~(e). 



B. Rectangular Billiard 

1. Spectral Staircase 

In Figs. |2^-|2j5, SS is shown respectively over a shorter 
and longer energy scales for several values of aspect ratio 
a. The former reveals noticeable oscillations around the 45° 
straight line. To further emphasize this point, in Fig. [Sfc, we 
plot jV(£) — £ for these a's. It is natural to anticipate that 
upon a-averaging, ^Y(e) — £ and the theoretical evaluation of 
5jV{e) using (Q} should vanish. However, numerical simula- 
tion shows neither to be the case, as seen from Fig. [2ji. While 
a translation of the latter downward and rightward bring the 
two into congruence, as seen from Fig. [2^, we do not fully 
understand the nature of this phenomenon IU6II . Clearly, it is 
inherent to the nature of parametric averaging and of the PO 
theory and is not an artifact of the numerical calculation. We 
remark that the proximity of (,jV(e)) — S to can be seen as 
a measure of performance of parametric averaging in attain- 
ing ensemble averaging. Comparing its magnitude to that of 
<y^(e) — e indicates that it does quite well. 



2. Global Level Number Variance and Correlation Function of 
Spectral Staircase 

In Fig. [3^, we plot E ? (g). Theoretical fit of the numerical 
data is quite good - given the limitations discussed above - 
and underscores importance of the non-diagonal terms. 



In Fig. [5J), we plot K^(e, co). Theoretical and numerical 
curves are in excellent agreement, which proves applicability 
of DA (0). We specifically point out the small CO behavior 
in the insert of Fig. [3K which is described by the small CO 
expansion in (0 and corresponds to the d(co) term in the level 
density correlation function. |Q1 



IV. SUMMARY 

We examined the global level number variance and the cor- 
relation function of spectral staircase in generic integrable 
systems with no extra degeneracies. We demonstrated that 
the global level number variance exhibits persistent oscilla- 
tions around the saturation spectral rigidity. These oscillations 
cannot be explained in the diagonal approximation framework 
and require an account of interference terms. 

Conversely, the correlation function of spectral staircase 
is well explained by the diagonal approximation. The latter 
points to the subtlety of the co — > limit since mathematically 
interference is destroyed by a finite co. 

In the future, we need to gain greater insight into integrable 
systems beyond the better understood rectangular billiards. 
For instance, we need to develop a quantitative description 
of the interference effects leading to the oscillations of the 
global level number variance in the modified Kepler problem. 
Larger-scale oscillations of the saturation spectral rigidity in 
circular and elliptic billiards, which also appear to be a prod- 
uct of the periodic orbit interference and are thus outside of the 



range of applicability of the diagonal approximation frame- 
work, are of great interest and call for further investigation. 
Properties of the spectral staircase, as well as its description 
using parametric averaging, require a closer examination as 
well. 
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FIG. 2: Rectangular Billiard: (a-b) Spectral staircase for six aspect ratio a's. Different colors encode different aspect ratios, (c) 

•j¥(e) — £ for six a's. (d-e) Blue line: JV (£) — £ calculated by averaging over 10 5 a's. Red line: theoretical 8JV{e) calculated 

from (Q]i and averaged over a's and in (e) the theoretical line is shifted 1 downward and 30 rightward. 





(a) 

FIG. 3: Rectangular billiard: (a) Comparison between numerical and theoretical L g , Black line: saturation spectral rigidity. 
Blue line: numerical £ g calculated from {[jV — e] 2 ). Green line: numerical L g calculated from ([jV — (.yK)] 2 ). Purple line: 

theoretical £„ calculated from (O after averaging over aspect ratios. Red line: theoretical £„ calculated from diagonal 

approximation plus interference between terms {M\ ,M-i) and {Mi,M\ ) with M\ ^ M%. (b) Correlation function of spectral 

staircase with e = 10 5 ; insert shows small 0) behavior. Blue line: calculated from {iyV{e\) — £\)(<jV(£2) — £2))- Green line: 

calculated from ((,_A^(£\) — (jY (fii)}) {^V{£q) — (,jV{£2))))- Purple: theory with diagonal approximation and parametric 

averaging. 



